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Abstract. Under certain conditions on k we calculate the limit distri- 
bution of eigenvalue number k, Xk, of the Gaussian unitary ensemble 
(GUE). More specifically, if n is the dimension of a random matrix from 
the GUE and k is such that both k and n — k tend to infinity as n — > oo, 
then Xk is normally distributed in the limit. 

We also consider the joint limit distribution of (x^^ , . . . , Xk„-^ ) where 
we require that ki, n — km and fci+i — ki, 1 < i < m — 1, tend to infinity 
with n. The result is an m-dimensional normal distribution. 



1. Introduction and formulation of results 

The Gaussian unitary ensemble (GUE) is a classical random matrix en- 
semble. It is defined by the probability distribution on the space of nxn 
Hermitian matrices given by 

F{dH) = Cne-^'''"'^^dH. 

By dH we mean the Lebesgue measure on the essentially different mem- 
bers of the matrix, namely 

(1.1) {ReHij; I <i < j < n,ImHij; I <i < j <n}. 

In other words this means that the entries in (|l.ip are independent Gaussian 
random variables with zero mean and variance i^^^. The measure on the 
matrices naturally induces a measure on the corresponding n real eigenvalues 
Xi. This induced measure can be explicitly calculated and its density is given 
by 



Pn{xi, . . . , Xn) — .g-) W 



,.2 ^2 



(2) Xji e 

l<i<j<n 

(2) 

The normalization constant Zn is called the partition function. It is often 
convenient to work with the eigenvalues being ordered. Naming the eigenval- 
ues so that xi < . . . < Xn, gives that the probability density pn,nixi, ■ ■ ■ , x„) 
of the ordered eigenvalues defined on the space 

^ord ~ {^1) • • • ) ^n] Xi < . . . < Xn} 

is given by 

Pn,n{xi, • • • ) ^n) — Pni^l, • • • ; ^n)' 
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This density p is a member of a family of functions called the correlation 
functions. These functions are defined by 

n\ f 

Pn.kixi, - ■ ■ ,Xk) = '-j^ / Pn{xi,...,Xn)dXk+l...dXn 

{n-ky.J^n-k 

= det{Kn{xi,Xj))ij^^. 

Here Knix,y) is given by 

n-l 

i^„(x,y) = ^/i,(x)/i,(y)e-5(-W) 

i=0 

where {/ii}i>o are the orthonormalized Hermite polynomials, that is 

2 

hi{x)hj[x)e~^ dx = 5ij. 



The kernel Kn{x,y) can also be represented by the so called Christoffel- 
Darboux identity. For a; 7^ y it holds that 

K (x v) = (^Y^^ hn{x)hn-l{y) - hniy)hn-l{x) ^_l^^2^y2^ 

" ' V2/ X — y 

and on the diagonal one has 

Kn{x,x) = (^nhl{x) - ^/r4JlTT}hn-l{x)hn+lix)J e"^^ 

The correlation function pn^i describes the overall density of the eigenvalues 
and the Wigner semi-circle law states that 

(1.2 hm A/-p„,i(V2nx = < ^ -fiiZ, 

n-+oo V n [0 it |x| > 1. 

All the results above and more can be found in the book by Mehta, [7]. 

This paper deals with the distribution of eigenvalue number k, Xk, of the 
GUE. More specifically, we look at the distribution of Xk as n and k both 
tend to infinity. For example, if 

k = k{n) = n — logn 

then, as n becomes large, k is very close (relatively) to the right edge of the 
spectrum. Another example is when k = n/2. In this case we are in the 
middle of the bulk of the spectrum. In both cases one ends up with a normal 
distribution in the limit. The following theorems generalize and specify this 
statement. 

Theorem 1.1 (The bulk). Set 

G{t) = - I - dx -l<t<l 
^ J-i 



GAUSSIAN FLUCTUATIONS OF EIGENVALUES IN THE GUE 



3 



and t = t{k,n) = G~^{k/n) where k = k{n) is such that k/n — > a G (0,1) 
as n ^ oo. If Xk denotes eigenvalue number k in the GUE it holds that, as 
n — > oo, 



^ ■ N(0, 1) 



1/2 



2)n 

in distribution. 

Theorem 1.2 (The edge). Let k be such that k ^ oo but - — > as n ^ oo 
and let Xn-k denote eigenvalue number n — k in the GUE. Then it holds 
that, as n ^ oo. 



J_^2/3 logfc \^/^ 



N(0,1) 



in distribution. 

Remark 1. The theorems deal with the bulk and the right spectrum edge. 
One gets the equivalent for the left edge with some obvious modifications. 



Remark 2. In [12] the distribution of the largest eigenvalue was studied. 
Remark 3. Set / = (— oo, s]. In [Ij it is shown that 



(1.3) P[xfc G (s,s + ds)] 
1 



Jfe(xi,...,Xfc_i,s);u(dxi)...;u(dxfc_i) /i(ds). 



(A; - 1)! Jjk 
Here is the so called Janossy density and 

//(dx) = C e 2 dx. 

In [1] it is also proven that Jk can be expressed explicitly by a determinantal 
formula. For k and n as in (II. ip or (jl.2p we should thus have that (jl.3p is, 
for large n, approximately equal to the probability density function for the 
normal distribution N{v,a). The parameters v and a should of course be 
taken to be those indicated from the relevant theorem above. 

Remark 4- The zero number k of the Hermite polynomial of degree n is 
close to the expected value of eigenvalue number k of GUE„. This can be 
shown directly by the following result, [4j: 
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There are constants ko and C such that ior ko < k < n — ko and a = k/n 



it holds that 

Zk,n 



(1.4) 



2n 



G- 



2-Kn 



arcsm 



{G-\k/n)) + 



2n 



< 



C 



n2(a(l -a))4/3- 

Here zi^n < ■ ■ ■ < z^^n are the zeros of the Hermite polynomial of degree n. 
When we are in the bulk this translates into 



Zk,n-V2nG ^{k/n) 



< 



C 



'n 



This means that one can replace t\phi by Zfc.n in Theorem ll.il Close to the 
edge, for example when k = n — logn, we can not use ()1.4p to conclude that 
this replacement is allowed. 

A motivation for this approximate equality between the locations of the 
zeros and eigenvalues goes as follows. Set 



W 



1 

1=1 



logi^i - 

l<i<j<n 

and note that 

Pn,nixi, ... ,Xn) = Const • e"^^. 

It is a fact, see that W obtains its minimum exactly when Xj = Zi^n, 
1 < i < n. This configuration is hence the most 'probable' for the eigenval- 
ues. Expanding around this minimum we see that it is reasonable that 
should have Gaussian fluctuations around Zk^n- 

Remark 5. If one is interested in the distribution of the eigenvalues of 
some other ensemble one should in many cases be able to apply the same 
methodology that we use in this paper. 

It is also interesting to see what happens when looking at several eigen- 
values at the same time. When we write k{n) ~ below we mean that 
k{n) = h{n)n^ where /i is a function such that, for all e > 0, 

h{n) 



(1.5) 



and h{n)n'' 



oo 



as n — > cxD. We have the following results: 

Theorem 1.3 (The bulk). Let {xki}^i be eigenvalues of the CUE such 



that < ki — fcj-i-i ~ n % < < 1, and ki/n 
Define Si = Si{ki,n) = G~^{ki/n) and set 



ai £ (0, 1) as n 



oo. 
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Then as, n — > oo, 

where is the cdj^for the m- dimensional normal distribution with covari- 
ance matrix Aij = 1 — max{0fc| i < k < j < m} if i < j and Aj^j = 1. 

Theorem 1.4 (The edge). Let {xn-ki}iLi eigenvalues of the GUE such 
that ki ~ n'^ where < 7 < 1 and < ki^i — /cj ~ n^S < < 7. Set 



Then as, n — > 00, 

F[Xi < Xi, . . . ,Xm < Xm] > ^a{xi, ■ ■ ■ , Xm) 

where ^\ is the cdf for the m-dimensional normal distribution with covari- 
ance matrix Aij = 1 — ^ max{^fc| i < k < j < m} if i < j and Ai^i = 1. 

Remark 1. As one would expect the eigenvalues near the edge of the 
spectrum are less correlated than the ones in the bulk. 

Remark 2. The eigenvalues are quite correlated in the bulk. In order for 
Xk and Xm to be independent in the limit it must hold that \k — m\ ^ n. It 
is interesting to compare with the following result by Mostellei0, [3] p. 201: 
Let Xi, z = 1, . . . , n, be independent random variables uniformly distributed 
on (0, 1). Consider the asymptotic joint distribution of the m sample quan- 
tiles Xn^ , j = I, . . . ,m, where Uj = [Xjn] + 1 and < Ai < . . . < Am- 

Theorem 1.5 (Hosteller). As n ^ 00 the joint distribution ofXn^ , ■ ■ ■ , ^rim 
tends to the m-dimensional normal distribution with means Xj, variances 
n~^Xj{l — Xj) and correlations 



Hence, in this case X^, ■ ■ ■ ,Xn^ are globally correlated in the limit. 



The proofs of Theorems 11.11 and 11.21 rely on a theorem due to Costin, 
Lebowitz and Soshnikov, see ^ and [iO\- Before presenting it we need some 
notation. 

Let {Pt}, t £ M+, be a family of random point fields, see [9j, on the 
real line such that their correlation functions have a determinantal fornji. 

^Cumulative distribution function 

Mosteller actually allowed for Xi to come from more general distributions. 
"^An example is the GUE. 





J < J • 



2. Proofs of Theorems 11.11 and 11.21 
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Call the determinant kernels Kt{x,y) and let {It} be a set of intervals. We 
denote by At the integral operator on L'^{It) defined by the kernel Kt{x,y), 
At : L'^ih) L'^ih)- Furthermore, Kt and Var^ denote the expectation and 
variance with respect to the probability measure Pt- Finally, let stand 
for the number of particles in It. 

Theorem 2.1 (Costin-Lebowitz, Soshnikov). Let At = Kt ■ Xh ^ family 
of trace class operators associated with the determinantal random point fields 
{Pt} such that \siTt{ifIt) = Trace(j4f — Af) goes to infinity as t ^ oo. Then 

\/Var(#/,) 

in distribution with respect to the random point field Pt . 

The following lemmas will be proven in sections H] and [5) 
Lemma 2.1. Let t = t{k,n) be the solution to the equation 



2 /■* 

n— / \/l — x^dx = k. 

71" J-l 



where k = k{n) is such that k/n — > a G (0,1) as n ^ oo. The expected 
number of eigenvalues in the interval 



In 



I — /log n \ 

V2nt + xy^^,oo 1 



is given by 



E[#/„] =n-fc--V(l-t2)logn + 



vr \ n 



log n 



Lemma 2.2. The expected number of eigenvalues in the interval In = 
[\/2n t, oo), where t — > 1~ as n ^ oo, is given by 

IE[M = ^n(l-t)3/2 + 0(l). 

Lemma 2.3. Let 5 > and suppose that t, which may depend on n, is such 
that —1 + 6 < t < 1 and n(l — t)^/^ ^ oo as n ^ oo. Then the variance of 
the number of eigenvalues in the interval In = [t^/2n, oo) is given by 

Var(#/„) = ^ log[n(l - t)3/2] (1 + n{n)) 

where lim„^oo vi^) = 0- 

Using these lemmas and Theorem l2.1l we are now ready to prove Theorems 
OandO 



Proof of Theorem Set 



wo , e I ^^S"- V''^ 
ty Zn + E — ;r- — , oo 

^ ' 4(1 - t^)n J ' 
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7 



Xk — t\/2n 



log n 
_ \^4(l-t^)n 



1/S 



log n 



4(1 - t^)n 



1/2 



\ [#/n <n-k]=¥„ 



(Var„ (#/„)) (y^^^(^j^))i/2 



(Var„(#/„)) 



1/2 



< C + e(n) 



where e(n) ^ as n — > cxd. By the Costm-Lebowitz-Soshnikov theorem the 
conclusion follows. □ 

Proof of Theorem \1.2[ Let g{t) be the expected number of eigenvalues in 
the interval In = [tV2n, oo). We have that 



Xfi—k 



< tV2n 



mn < k] 



#ln - git) 



< 



k - git) 



If we can find t such that 
(2.1) 



.(Var„(#/„))i/2 - (Var„(#/„))i/2 
k - git) 



(Var„(#/„))i/2 

as n ^ oo then, by the Costin-Lebowitz-Soshnikov theorem, we are done. 
The idea now is therefore to find a candidate for t. We will then insert this 
t in the equation above to see if it is satisfied. Set for simplicity /i(t) = 
(Var„(#/n))i/2- We get from Lemma [2^ and Lemma [2^ that 

g(t) =ain(l-t)3/2 + 0(l) 

hit) = 02 logl/2[^(i _ ^)3/2] ^ o(logV2[„(l _ t)3/2]) 

where Oj are known constants. We need to study the equation 

k = git)+Chit) 
or, since 5 is a strictly decreasing function, 

t = g~\k - Chit)) ^ g-\k) - ig-'Yik) ■ ^hit). 

Since 

ig-'Yik) 

we need to study g^^ik). 



1 



g'ig-'ik)) 



k « ain(l - 

^ X 2/3 



1 



ain 
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A reasonable guess for the derivative of g is that 



9 it) 



3a 



-nVl-t. 



We now get 



1/2 



and 



h{t)^h{g-\k))^a2log^/^ 
When gluing the pieces together one gets 

/ A- 
t-l-f — ) 
ain I 



n- 



.:^^l/3^2/3 



a2 log^/2 h 



2a2 log^/^fc 

2/3 /^l/3^2/3 • 



3a^ 

When inserting this expression in (j2.ip it turns out that it all works out. 
Some rearranging finally yields the result. □ 



3. Proof of theorems 11.31 and 11.41 
We shall use the following theorem, [llj: 

Theorem 3.1 (Soshnikov). Let {X,J^,Pl) be a family of determinantal 
random point fields with Hermitian locally trace class kernels Kl and 



Ul '•••i^L SL>0 

be a family of Borel subsets of M, disjoint for any fixed L, with compact 
closure. 

Suppose that the variance of the linear statistic Yli^-oo fL{xi), where 



fiix) = ^ • XjU) {x) ai, . . . , Ofc G M, 
grows to infinity with L in such a way that 



(3.1) 



Vari(#4^'))=o[varz.( /l( 



for any i < j < k. Then the Central limit theorem holds: 



7V(0,1) 



in distribution. 
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Remark 1: The theorem in is actuaUy more general than the theorem 
stated here. 

Remark 2: If (|3.ip holds for any Q!i,...,afc then , . . . are 

jointly normally distributed in the limit, p]. 

Proof of Theorem \1.3[ Take ki, Si, 9i and Xi as in the formulation of The- 
orem [T31 Note that ki — fcj+i ~ n^' implies that Sj — Sj+i ~ vP^~^. 
For any real numbers Xi we have the identity (for n large enough) 



^ #/i - E[#/i] ^ n - fci - E[#/i] 
[(Var(#/i))i/2 - (Var(#/i))i/2 ' 

#Ji + #/2 - E[#/i + #/2] ^ n - fc2 - E[#/i + #J2] 
(Var(#/i + #/2))V2 - (Var(#/i+#/2))V2 

(Var (Er=i #^.))'/' (Var (^"=1 



where 2 < i < m. We will now show that the random variables 



are jointly normal in the limit. To this end we shall use Theorem 13.11 above 
to show that all linear combinations of the variables are normally distributed 
in the limit. Since 



Here the intervals /, are given by 




m 



#/l,#/l + #/2,...,^#/i 



i=l 



Qi#/i + a2(#/i + #/2) = (ai + a2)#/i + a2#/2 
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and so forth it is clear that one can instead look at all linear combinations 
of Hence, by Theorem 13. lln we must calculate (see Appendix B) 

Var (ai#/i + a2#/2 + • • • + am#Im) 



m „ „ 

m „ „ 

- ^ aiUj I i K^ix, y) dxdy 



to see that it is of magnitude log n. First define the set M hj k £ M <^=^ 
0fc = 1. Hence 

M = {ki, . . . , kj}; l<ki<k2<...<kj<m-l 

for some j such that < j < m — 1. Suppose first that M is empty which 
means that Oi <1 for all i. If ai 7^ then, by using the inequality 

xy < ^(x^ + 

we get 

Var (ai#/i + a2#h + • • • + am#Im) 

>Y.'^U ^n(^,y)dxdy 

i^j JJliXlj 

(3.2) =E"?(// ^n(a^,y)dxdy 

-fill Kl{x,y)<l^dy\ 

All the terms under the z-summation in (j3.2p are non-negative and the first 
term can be calculated as in the proof of Lemma 12.31 which we provide in 
section [5j In the proof of this lemma we show that in the domain 

^ = I {x , y); s < X < s + , s < y < s 

[ log n log n 



it holds that 



2„if„( V27,y) = + " ( 



^The theorem by Soshnikov does not apply directly to this situation since Ji does not 
have compact closure. This is however easily overcome simply by chopping of the interval 
far out where the probability of finding any eigenvalue is exponentially small in n. 
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It is also shown that if 

fi' = y); V^2ns < X < oo, — cx) < y < \/2ris| / V^n ■ Q 

therS 

// Kl{x,y)dxdy = 0{\oglogn). 
J Jn' 

In what follows we shall often make use of these facts without mentioning 



it. The main contribution to the first term in (|3.2p can now be calculated 
to be (disregarding af) 

J 1 u _ y\2 ^y^^ = ^ + log n) 

^ log n 

where 9* = maxj Oi <1. By our definition of ~ above the integration in the 
y- variable should have been over the interval (si — 1/logn, si — h{n)n^ 
where h{n) satisfies (|1.5|) . However, because of the logarithmic answer this 
h will only produce lower order terms. 

Now suppose that j = as before, ai = . . . = ak-i = but Ok / 0. In 
this case we get 

Var {ak#Ik + . . . + am#Im) 

^fl'^H If Klix,y) dxdy 



f2 [[ Klix,y)dxdy), 



Using the estimates above it is straightforward to verify that the /c-term is 
of order logn. 

When J > 1 meaning that there is at least one k with 0^ = 1, things are 
only slightly more complicated. Let k* be the largest integer i such that 
6i = 1. It is sufficient to consider the case when there exists i > k* + 1 such 
that 7^ 0. On the other hand if this is the case then we are in a situation 
very similar to when j = 0. Either ak*+i / or ak*+i = . . . = a/_i = 
but ai 7^ 0. The details are left out. 

It is hence a fact that 

m 

#Il,#/l+#l2,...,^#Ii 

1=1 

in the limit have a joint normal distribution. 

To complete the proof we need to calculate the correlations between the 
different #/i's. If j < iwe have that Sj—Si ~ n~"' where 7 = 1— max,<fc<j 9k- 



^Since K„{x,y) = Kn{y,x) it is clear that the same estimates hold in the domains 
obtained from reflection with respect to the x = t/-line. 
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Set 

k 

m=l 

From a straightforward calculation (as above) we get that 



Var(X, - Xj) = Var I #4 = Var # J 4 

fc=j+i / \ k=j+i 



1-7 

— ^ log n + O (log log n) . 



Since 

Var(Xfc) = — log n + ©(log log n) 
the correlation p is given by 



i (Var(Xj ) + Var(X,) - Var(X j - X^-)) 
VVar(Xi)Var(Xj) 

□ 



Proof of Theorem \1.4\ This proof is of course very similar to the previous 
one so some details will be skipped. 

With notation as in the formulation of Theorem 11.41 the intervals of in- 
terest (confer with the previous proof) are in this case 

1/2 



\ 1/2 

log h-i \ 



where Ci, C2 are known constants and 2 < i < m. Given any xi, . . . , Xm it 
is easy to see that the sets Ij are intervals if n is sufficiently large. As in the 
previous proof we want to show that 

are jointly normally distributed. The way to prove this is the same as before 
but some details are different. By Lemma 12.31 we need to show that 



logn = Var K;;ai#/, 



for any real a^'s such that Qj / for some i. 
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Let t = t{n) be such that n'^~3 < 1 — t < for some < e < 1/3. From 
the proof of Lemma 12.31 below one can see that in the set 

Qt = { x,y ; t < X <t + - ,t-- <y<t 

I log n log n 

it holds that 

Returning to the variance calculation we first assume that ai ^ 0. We know 
from the previous proof that in this case it is sufficient to to show that 



JJhxd'rW I"!. I,) 



is of order logn. In fact since the integrand is non- negative it is enough if 

1 



■ dydx 



is of order log n where 



and 



r = {ti+ri,ti + \ 

log n 

U — [ ti — ,ti 

logn 



XiC2 



log h \ 



An elementary calculation shows that this integral is indeed of order logn. 
If ai = . . . = Ofc-i = but 7^ it is sufficient that the integral 

dydx 



is of order log n where 

J = [ tk-i + rk-i,tk-i ^ : 

\ log n 

J* = {tk, tk-i) ■ 

Again we get the size logn. This proves that we get a normal distribution 
in the limit. The calculations of the correlations are very similar to the bulk 
case and the details are not presented here. □ 
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4. The expected number of eigenvalues in I„ 



In this section and the next we shah need asymptotics for the Airy func- 
tion and the Hermite polynomials. In [4] the asymptotics for a class con- 
taining the Hermite case was studied. It is shown that there exists a 5 > 
such that the following holds: 
1. -l + 5<x<l-5 



2nx e 



^^2^ (1 - x2)l/4 



COS 



2nF{x) — — arcsin(x) 



2. 1-6 <x <1 



(^■\/2nx^( 



(2n)-i/^ 

1/4 



{ (r^) [3nF(x)]^/^ Ai (- [3nF{x)f') (1 + ©(n-^)) 



1 — X 
1 + x 

3. l<x<l + 5 



x + 1 



1/4 



[3nF(x)]'^/^ Ai' ( - [3nF{x)f/^ ) (1 + 0{n~^)) 



2/3 



hn{\'2nx)e 



{2ny 

1/4 



-1/4 



x-l 

X — 1 

x+1 



[3nF(x)]^/^ Ai([3nF{x)] 



l2/3 



1/4 



[3nF{x)r'^/^ Ai' ([3nF{x)f^) \ (1 + 0(n"i)) 



4. x> 1 + 5 



hn{v 2nx)e 



O ( n-i/4e-"^(-) 



In these expressions Ai stands for the Airy function and 
(4.1) F{x)= /Vll-y'ldy 

There are of course also similar asymptotics for the Hermite polynomials 
near the point —1. 

The Airy function is bounded on the real line. It is exponentially small 
in X on M_)_ and for r > it holds that, [8j, 



Ai(— r) = tt' 



"1/2^-1/4 



COS 



.3/2 _ ^ 



Ai'(-r) = vr-i/Vi/^ 



sm 



f 3/2 _ ^ 

3 4 



+ 0(r-3/2)| 
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Proof of Lemma [2A[ Set 



fnit) =t + X 



2n 



n 



We have that 



/n{t) 



npn{x) dx 



where pn is the scaled density for the eigenvalues (the limiting density has 
support in [—1, 1]). From symmetry one gets 



fn{t) 



npnix) dx 



n 



/n(t) 



npn{x) dx. 



Formula (4.2) in ^ applied to the Hermitian case says that for 5 > it 
holds that 



np 



'n(x) =n—\/l — 
vr 



if X e [-1 + 5, 1 - (5]. We get that 

„ 2 /'•/'"(*) , 

^[#In] = ^-n- Vl-x^ dx+0(n~i) 

^ Jo 

2 /■/"(*) 



+ 0(n 



n — n- 



rlnW 

J Vl-x2dx+0(n-l) 



n — n 



n-k- -yJ{l-t'^)\ogn + 

IT 



log n 



n 



□ 



Proof of Lemma \2. ^ From the formulas (4.4) and (4.21) in [6j one gets after 
some minor calculations that 



npn{x) 



<I>'(x) 7'(x) 



[2Ai($(x))Ai'($(x))] 



4<I)(x) 7(x) 
+ $'(x) [(Ai'($(x)))2 - $(x)(Ai($(x)))2] + O 



n 



in a fixed neighborhood of [0, 1]. As before pn is the scaled density for the 
eigenvalues so that 



npn{ 



x) dx 
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and the functions 7 and $ are given by 

IX 1/4 



X 



x + l 



- dy 



(Sn/i^V^dy) 



2/3 



2/3 



if X < 1 
if X > 1. 



The function 7 is evaluated by taking the Umit from the upper half plane 
using the principal branch. 

The fact that the asymptotics only holds for x G [0, 1 + (5] for some 5 > 
(independent of n) is not a problem. It is not difficult to show that, for 
x > 1 + 5, Pn{x) is exponentially small in n and exponentially decaying in 

X. 

We shall now study the different terms in the asymptotical expression for 
Pn above. When looking at the asymptotics for Ai and Ai' it easy to see 
that 

|Ai(x)Ai'(x)| =0(1). 

A calculation shows that 

$'(x) iixY 



hence 



4<I)(x) 7(x) 
^+^(^'{x) i{x) 



0{l) 



n 



[2Ai($(x))Ai'($(x))] dx = 0{1). 



The main contribution comes from the second term. In fact a primitive 
function can be found for this expression (see Appendix A): 

.1+5 



$'(x) [Ai'^($(x)) - $(x)Ai2($(x)) 



dx 



[set y = <l>(x)] 



Ai"(y)-yAi^(y)dy 



^(yW(y)-2/Ai'^y) 



-Ai(2/)Ai'(y) 



--Ai($(t))Ai'(<I>(t)) + 0(e— ) 

Here c is a positive constant. Integrating the third term only gives a contri- 
bution of order n~^. One can now use the asymptotics for the Airy function 
and its derivative to get the stated result. □ 
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5. The variance of the number of eigenvalues in /„ 



Proof of Lem'ma \2.S^ The proof will be divided into two basic cases. The 
first case is when l — t>6 for a fixed (5 > 0, that is, we are in the bulk. The 
second case is when t = t{n) — > 1~ as n — > oo, that is, we are close to the 
right spectrum edge. 

First define = [tV2n, oo) and #1^ as the number of eigenvalues in /„. 
It is a fact (see Appendix B) that 



Var(/„) 



K^{x,y)dxdy 



T / 



where Kn is the Hermite kernel. The advantage with this representation 
is that there is only one singular point involved in the Christoffel-Darboux 
representation of Kn{x,y): 



Kn{x,y) 



n hn{x)hn-l{y) - hn-l{x)hn{y) ^_l(^^2+y2) 



-e 2v 



x-y 



Case I (the bulk): After a change of variables (x \^2nx) we get the 
integrand 

\/2n Kn{V2nx,V2ny) 
First consider the domain where both variables are in the Bulk: 
(5.1) F = {{x,y);t< x<l-5,-l + 5<y<t} 

Recall that in F, /i„ has asymptotics given by 

/i„(^/2^x)e-"^' 



1/2 



1 



where 



7rv2n 



F{x) 



(l-x2)l/4 



COS 



2nF{x) — - arcsin(x) 



+ 0(n-i)^ 



j a/i — dz = - ^arccos x — x\l 1 — 



The asymptotics for hn-\ becomes 



1/2 



/i„_i(\/ 2nx)e 



(l-x2)l/4 

2 



7rV2(n-l)y 

2(n - \)F(xn) - ^ arcsin(x„) + OirT^)^ 



cos 



7r^/2^y (l-x2)i/4 

X ( COS 



2(n — \)F(xn) — — arcsin(x) 



+ 0{n 
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where Xn = \ '^^^x. A Taylor expansion gives 



F{Xn) = F{x) 



2(n- 1) 



- x2 + C)(n-2) 



and hence 



2(n - l)F{xn) = 2nF{x) - 2F{x) - x\/l - x2 + 0{n~^) 

= 2nF{x) — arccosx + Oin^"^) 

One can now write 
hn{V2^x)K-i{V2^y)e-<'''+y'^ 
2 



7r\/2^(l-x2)l/4(l_y2)l/4 



X COS 



2nF{x) — — arcsinx 



cos 



2nF{y) — - arcsin y — arccos y 



Set 



ax = 2nF{x) — — arcsin x 



6x = arccos x. 
By the Christoffel-Darboux formula 



2n Kn{V2nx, \/2ny) 



1 



7r(l-x2)l/4(l_y2)l/4 

cos ax cos \ay — Qy[ — cos [a^ — Qx\ cos + O (n~ -"^ ] 
x-y 

To prepare for integration we now divide F into four disjoint sets. Set 







< X < t + -,t <y <t 

n n 



Ti = T\\JTl = { {x,y);t < X <t+^,t- ^ <y< 
' r[n) r[n) 



(x,y);i + - < X <t + —-,t- - <y<t 
n r[n) n 

= r\(rouri) 

where r(n) = logn and T was defined in ()5.ip . 

Fq: When integrating over Fq one can use the fact that 

sin(x — y) 



'2n \KniV2nx, V2ny)\ < Cn 
where C > 0. Hence 



x-y 



2n Kn{V2nx, V2ny) dxdy = 0(1 
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Fi: For x G Fi it holds that 

9x = arccos x = arccos t + O 



1 



r(n) 



and of course also the equivalent for 6y. Defining 9 = arccos t we get by the 
use of some trigonometric identities that 

cos ax cos[ay — 9y] — cos[aa; — 9x] cos ay 

= cos ax cos[ay — 9] — cos[q;j; — 9] cos ay + O ( — ) 

\r{n) / 

= Vl - sm[ay -ax] + o( ) . 

\r{n)J 



Since 



and 



(l-x2)l/4(l_y2)l/4 



r(n) 



ax-ay = 2n{F{x) - F{y)) + O (^-^-^ 



we now get 



2nK„(v2nx, V 2ny) dxdy 

1 sin2[2n(F(x)-F(y))]+0(^ 



^ (x — uY 

l-cos[4n(F(x)-F(y))] 



7!"' 



dxdy + 



//. 



{x - y)2 

cos [4n{F{x) - F{y))] 



dxdy + O (log r(n)) 



dxdy 



{x - yf 



dxdy + O (log r(n)). 



The remaining integral is not bigger than a constant as will now be shown. 
A partial integration in the y-variable gives 



cos [4n(F(x) - F(y))] 
{x - y)2 



dxdy 



sin [4n(F(x) - F{y))] 
AnF'{y){x - y)'^ 



.1+1 
r(n) 



sin [4n(F(x) - F{y))] 



4n[F'{y){x-yr]Jy 



dy dx 



=: h - h. 
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Both these integrals are easy to estimate: 

1 



We have 



n{x — {t — n 

y 



dx = 0(l). 



[F'{y)ix-yf] ^) - , 

' Jy {l-y^yi\x-yY ^i-y2(^x-yY 

which gives 



I/2I < C 



rr _^ 

JJr\ n{x - 



yf 



o{i). 



T2: In T2 it holds that 



2nKn{V2nx, ^/2ny) 



O 



{x - y)2 



and a trivial calculation gives 



1 



■ dxdy = O (log r(n)) . 



'J 1-2 (-^ - yy 

To complete case I we must also integrate over x \ F. By using the 
appropriate asymptotics for the Hermite polynomials, the Airy function and 
its derivative and then taking absolute values in the integrals it is straight 
forward to show that the contribution from this domain is 0(1). We omit 
the details. 

Case II (the spectrum edge): First consider the subdomain 

Vt= {{x,y);t<x <l- Cn~^, 1-S<y<t} 

where C is a large positive constant. After a change of variables the contri- 
bution Jci from \/2n • Q to the variance can be written as 



Jn 



IIci 



2nKn{y2nx,v2ny) dxdy 



In order to deal with this integral we must first study the integrand and, 
via Christoffel-Darboux, especially the difference 



(5.2) D 



hn V 2nx hn 



2ny 



1 (^\/2nx^ hn (y2ny^ 
We will show that in 0, it holds that 

b{4:n{n - l)y/^D = Ai [3nF{x)f^^ Ai' [3nF{y)f^^ 

-Ai' [3nF{x)f^^ Ai (- [3nF{y)f^^ 



+ 



1 



n(l — x) 



+ 



(1 - X)V4 
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Here 6 is a constant, Ai stands for the Airy function and 

F{x) = j Vl - t2 dt. 

J X 

In n hn has the following asymptotics: 

1/4 

I 

1/4 



(2n)-i/^| [3nF(x)]i/'^Ai (-[3nF(x)]2/3) (1 + ©(n-^)) 

\ 1/4 ^ 

^ j [3nF(x)]~i/6Ai' (-[3nF(x)]2/3) (1 + ©(n-^)) I 



1 — X 



l + x 



If we disregard the 0{n ^) terms in the /i^— asymtotics, then ()5.2p can be 
written as a sum of four differences Di — 



l + x\ / l + yn\ ro^' p^, ^^ Ve 



(4n(n - l))i/^Z)i 

/I , ^\ 1/4 /I , \ 1/4 

, . , [3nF(x)]^/^ [3n'F(y„)] 

X Ai (- [3nF(x)]2/3) Ai (- [3n'F(y„)] 

--I I \ 1/4 /I 1 \ 1/4 

'^'^"^ ^[^) [3n'F(x„)]^/^[3nF(y)]i/« 

X Ai (- [3n'F(x„)]'/') Ai (- [?,nF{y)f''' 

(4n(n - l)f'^D2 

X Ai (- [3nF(x)]2/3) Ai' (- [3n'F(y„)] '^^^ 
X Ai (- [?>n'F{xn)f'^) Ai' (- [3nF(y)]2/3 

(4n(n- l))^/^i:'3 

1 \ 1/4 /I I \ 1/4 

'-""^ + [3n'F(.„)]-^/^[3nF(y)]i/^ 



X Ai' (- [3n'F(x„)]'/') Ai (- [?,nF{y)f^) + 
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X Ai' (- [3nF{x)f^) Ai (- [3n'F(y„)]'/') 



(4n(n - l))^/^i^4 

X Ai' (- [3nF{x)f') Ai' (- [3n'F(y„)] 
X Ai' (- [3n'F{xn)f^) Ai' (- [3nF(y)]2/=^) 



We have used the notation n' = n — 1 and Xn = y:^^x. Note that Xn < 1 
in 17. 

Di : A calculation using the series expansion 

^ ' =co + ci(l-a;) + ... 



(1 - x)V4 
gives 



= (t^)''' (r^)''' t'"''^^^]''' 



where 



Since 



ai = hm Vl + x-^^ ^ . 



^'(-'='"^W'''')=°( „i/a(l-.)V^ 
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it holds that 

(4n(n - l))^/^L>i 

,1/3 



ain 



Ai (- [3nF{x)f/') Ai (- [3n'F{yn)f') 

-Ai (- [3n'F{xn)f') Ai (- [3nF{y)f^) 

■(l_y)3/4- 



+ 



(1 - a;)V4 y ■ 



D2 — D4: The same procedure as above gives 
(4n(n- l))i/^L>2 



Ai (- [3nFix)f^) Ai' (- [3n'F(y„)] 
-Ai (- [3n'F{xn)f') Ai' (- [3nF(y)]2/=^) 

'(l_y)5^- 
(1-X)V4 



+ 



(4n(n- l))i/^L»3 
= 0(1) 



Ai' (- [3n'F{xn)f') Ai (- [3nF{y)f^) 
Ai' (- [3nF(x)]2/3) Ai (- [3n'F(y„)]'/' 



+ 



■(l_y)5^^ 
(1 - X)V4 



(4n(n - 1))^/^L»4 



0{n 



-1/3^, 



Ai' (- [3nF(a;)]2/3) Ai' (- [3n'F(y„)] 
-Ai' (- [3n'F{xn)f') Ai' (- [3nF(y)]2/3) 

Now consider the difference still left in Di: 

Ai (- [3nF{x)f^) Ai (- [3n'F{y^)f^) 

- Ai (- [3n'F{xn)f') Ai (- [3nF{y)f' 
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To deal with this expression we shall first investigate the argument 

2/3 

,2/3 / , X ' 



3(n- 1)F 



n — 1 



[3n'F(x„)]- 
A simple integration shows that 

F{x) = J ^/l — t^dt = ^ ^arccosx — x^/Y^^x^^ . 



Since 



we get 



n 



n-1 



X = X + 



2{n-l) 



F{xn) = Fix) + F'{x) + + 0{F"{x)n- 



= n^) - 



xVl — x'^ 
2{n - 1) 



n- 



and hence 



Sn'F{xn) = 3(n - l)-F(x) - - + O I — ) 

2 \?^v 1 — X J 



= 3nF(a;) — - arccos x + 



ny/1 — x 



The argument can now finally be written as 
(5.3) -[3n'F{xn)] 



■ + 



(3nF(x))V3 Vn^/^(l-a; 



Note in the last expression that 

arccos x 



n 



-1/3 



(3nF(a;))V3 

It is now possible to expand the difference in a Taylor series around the 
point — [3nF(x)]^^'^ and the result is 

(12 



n 



1/3 



Ai [3nF{x)f^^ Ai' (- [3nF{y)f^^ 

- Ai' (- [3nF(x)]^/^) Ai (- [3nF(y)]2/''^ 



where 02 is defined by 



+ 



02 = lim 



(1 - y)'/" 

n4/3(l-x)y ' " \ nV3(i_a;)iAi 



+ 



arccos x 



X-.1- (3F(a;))V3- 
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Similar computations can be done in D2 - -D4 and one then ends up witli 

(4n(n - 1))'/\D2 + D3 + D^) = O ((1 - y)^^) . 
Adding everything up we now have that 
(4n(n - 1))V4 



-D 



aia2 



(5.4) 



Ai [3nF{x)f^^ Ai' [3nF{y)f'' 
-Ai' (- [SnF{x)f^^ Ai (- [3nF{y)f^^ 



+ 



n(l — 



+ 



(1 - X)V4 



As we shall see the main contribution comes from the domain 

ill = s {x,y);t < X < t -\ TT' * < 2/ < ^ — e 



r(n) ' r(n) 

Here r(n) is a function tending slowly to infinity as n tends to infinity and 
e = J the size of the expected distance between two eigenvalues at t. 

The reason why this e is necessary lies in the asymptotics for the Hermite 
polynomials. The error term given there, however small, will cause problems 
since the integral 

is divergent. 

From the asymptotics of the Airy function and its derivative we get, for 
(x, y) G r^i, that 

Ai [3nF{x)f/^^ Ai' [3nF{y)f/^^ 

= (^-^{nF{x))-^/^ sin ^2nF{x) + j\+0 ((nF(x))-^/6)^ 

X (^-^(nF(y))i/6sin [2nF{y) - ^] + O (^{nF{y)y'/')^ 

2nF{x) + ^1 sin \2nF{jj) - ^1 + O ((nF(x))"^) . 



1 (F{y) 



1/6 

If we define r(n) by 



sm 



r(n) 



= max 



(^^^'logKl-t)3/2]) 
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Fix)) 

1 

n(l — x) 

(l_y)3/4 



l + O 
l + O 

o 
o 



1 



r(n) 
1 

r(n) 



1 



nF{x 
1 



r(n) 



r(n) 



(l-x)V4 

Prom this it follows that in D can be written as 

(4n(n - 1))V4 



-D 



aia2 



-—[ sm 

TT 



— sm 



2nF{x) 


vri 
+ 4. 


sin 


2nFiy) - 


TT' 

4. 


2nF{x) 


TT" 
~ 4. 


sin 


2nF{y) + 


TT" 

4. 



=i sin [2n{F{x) - F{y))] + O f-^ 



The numerator in the integral of interest is 



n „9 



2^n{n - 1) 



+ 



4^2 ■sin^[2n(F(x)-F(y))] + o(-l^ 

= _ cos [4n(F(x) - F{y))\) + O 



We used here that 0102 = 2. A simple integration gives 

/ / , ^ ,2 = log[n(l - tf] + 0(log r(n)). 

The integral 



f2i 



COS [4n(F(x) - F{y))] 
[x - yY 



dx dy 
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is By doing a partial integration / can be split into two integrals: 



t-€ 



sin [4n(F(x) - F{y))] 
—AnF'{y){x-yY 



sin [An{F{x) - F{y))\ 



t-e 



t- 



r(n) 



t- 



r(n) 



AnF'{y){x-yyjy 



dy dx 



- + /2 



We can estimate 7i in the following way. 



iAi<2y^ 



t 4nVl - - (t - e))- 



dx 

e 
' 2 



X — t + e 



t+- 



{") e 2 



Since 



([ny)(^-y)r')|^ 



we get 



\h\<c(ff 



n(l - y)3/2(a; - y)2 ^^^"^ //^^ n-y/l - y{x - yf ^^^^-^ J 



■ dxdy ) . 



The first part is small: 



y)3/2(a; - y)2 
1 

< 



1 



The second part is also easily estimated: 



dxdy = O 



'log [n(l - tfl^] 
n(l - 1)3/2 



1 nVl - - y)- 



dxdy<e / /^^^— -^dxdy = 0(l) 



This concludes the calculations in f^i. 

The calculations made above can also be applied to the thin slice 

1 - t 

{x,y);t + e<x <t-\ -^,t — e<y<t 

r{n) 

and the result is O(log[r(ri)]). 
The corner 



r^o = {{x, y)',t < X < t + e,t — e < y < t} 
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requires a special technique. In this domain a different representation of 
will be used, namely 

i=0 

By use of the Cauchy-Schwartz inequality we get that 

Kl{x,y) < Kn{x,x)Kn{y,y). 

Having separated the variables we can now use the calculations from the 
proof of Lemma 12.21 to see that 

^ (y2^Kn{V2^x, V2^y)y dxdy = 0(1). 

Note that 



rt+e 

J Kn{V2nx,V2nx)dii = g{t) - g{t + e) 

where g{t) is the expected number of eigenvalues in the interval (t\/2n, oo). 

Now we shall look at the other part still left of 0,. This domain can 
conveniently be written as 0,2 U f^s where 

n.) = I {x,y);t < x < 1 - Cn~^, 1-5 <y<t 



and 



r(n) 



(x, y); t + -7-^ < x < 1 - Cn'^,t - -j-^ <y <t 



r{n) ' r{n) 

When looking at the expression for D in (|5.4p above it is clear that every 



term is smaller than 
n-i/^Ai [_ [ZnF{x)f^) Ai' {- [ZnF{y)f^ 



O n 



-1/2 



\ — X 



l/4^ 



This means that it is sufficient to calculate the integrals 

— =^^^-— dxdy i = 2,3. 
Vl — x{x — yj^ 

The calculations are straightforward so some details will be skipped. When 
first integrating with respect to the x-variable one gets 



1 VI - x{x - y)' 



■dx 



2(1 -y) 



log 



(yr^ + yi^) ( - yr^) 
1/1 1 
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Q.2'- Letting Hi = 1 instead of 1 — Cn^^ wc get nicer expressions. This is 
allowed since the domain of integration becomes larger. The task is to get 
an upper bound for the integrals 



A 



rH2 

JLn 



log 



dy 



-log 



and 
B 



dy 



dz 



2 r ' ( I I 1 

JVl^ \Z- Li Z+y/1- LiJ 



dz 



where 



L2 = l-S, H2 = t- 



r{n) 



and Li = t . 



When manipulating the integrand in A one gets 



log 



1 + 2 



1 



-O 



Some algebra shows that 



z{z- VT^) 
which can easily be integrated: 

'z-V^-Li 



A<C 



log 



The integral B is even easier and one gets 



B 



log 



Z-y/T^l 
Z+y/l-Li 



O [log r{n)) 



O {log r{n)). 



fl^: The same procedure as in gives that the contribution to the variance 
from this domain is o(l). 

We shall now consider the thin strip 

^4 = {x, y;l- Cn~'^ < x < 1 + Cra~\ 1-S<y<t}. 

The asymptotics here is similar to that in 0, and hence many of the calcu- 
lations already done can be applied here as well. As before D can be split 
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up in Di — D4 which can all be treated similarly. Therefore we only look at 
Di here. We have that 

(4n(n - l)f'^Di = ain^^'^ Ai [2,nF{x)f'^'^ Ai [3n'F(y„)] 

-Ai U [3n'F(x„)]'/') Ai (- [?,nF{y)f^ 



+ 



(l-x)l/4 1 ■ 



where 

T [3nF(x)] 



|2/3_/ - [3nF(x)]2/=^ ifx<l 
[3nF(x)]^/^ ifx>l. 



This follows from the calculations done above and the asymptotics for the 
Hermite polynomials when x > 1. In 04 we have 



Ai [3nF(x)]2/3j = Ai(0) + 
Ai(^[3n'F(x„)]'/') = Ai(0) + O(n-i/3) 
and by using equation ()5.3p (for the y-variable) one gets 

,4„(„-l))V^«, =0(115!;;), 



The error term here has actually already been dealt with in the estimations 
of the contribution coming from 

Rather than to repeat a lot of calculations we now just give ideas of how 
to to treat what is left of [t,cxD) x (— oo,t]. 

In the domain 

{x, y; 1 + Cn'^ < x < I + 5,1 - 5 < y <t] 

one can perform much the same calculations as in Vt and the contribution is 
In 

{x,y]t<x<l + 5,-l-5<y<l-5} 

one can use the fact that x — y > 5 to show that the contribution from this 
domain is 0(1). Ifx>l + 5ory<— 1 — 5t one easily gets from the asymp- 
totics for the Hermite polynomials that ivr„(\/2nx, \/2ny) is exponentially 
small in n and exponentially decaying in (or y^). Thus the contribution 
from this domain is o(l). □ 
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APPENDIX 

A: Some integrals. The following equalities hold: 



l-OO 

/ M\y)dy = 

Jx 








poo 

/ yAi2(y)dy = 

Jx 


- (xAi' (x) - a;^Ai^(x) - 


Ai(x)Ai'(x)) 




/ Ai'(y)dy = 

Jx 


^ (x'^Ai\x)-xAf{x)- 


2Ai(x)Ai'(x)) 




POO 

Jx 


^ (x^Afix)-x^Ai\x)- 


- 2a;Ai(x)Ai'(x) 


+ Ai2(x)) 


POD 

/ yAi''(y)dy = 

Jx 


^ (^x^Ai\x) -x'^Af{x) 


- 3xAi(x)Ai'(x; 


1 + ^Ai2(.)) 


The first integral is 
membering that 


obtained by performing a 


partial integration while re- 



Ai (x) = .TAi(x). 

The integrals 3-5 can be obtained rather easily from the second which can 

be treated as follows: 

Set 

Ua{x) = Ai(Q:x) a > 0. 

The relationship 
holds since 

u'^{x) = a^Ai"(ax) = a^xAi{ax) = a^xua{x). 

Hence 

POO 

/ XUaix)Ui3{x) dx = 3 _ 3 [u'^Up - Uau'p] ^ 
J a Oi p 

_ Ua{a)u'p{a) - u'^{a)up{a) 

The idea now is to let a, /? tend to one. Set a = 1 + h and (3 = l — h where 
/t > and small. The left hand side tends to 



POO 

/ a;Ai^(a;)dx 

J a 



as /i ^ 0+. Standard calculations show that at the same time the right 
hand side tends to 

^ [-a^Ai^{a) - Ai(a)Ai'(a) + aAi''(a)) . 
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B: Variance calculations. Let Ii, . . . , be a set of disjoint intervals and 
#/i be the number of eigenvalues of the GUE„ in the interval /j. We shall 
give a formula for Var {ai#Ii + . . . + am#Im)- Clearly 



n 



k=l 



where xb is the characteristic function for the set B and {xk}k=i ^^re the 
unordered eigenvalues. The expected value is easy to compute: 



E 



Jl J h 



The correlation functions Pn,k were defined in the introduction. We also 
need to calculate E[#/?]: 



E 



j,k=i 



^nxhixk)]+Y,nxhixk)xh{xj)] 



k=l 



' Kn{x,x)dx+ // p„,2(a;,y)dxdy 

L JJliXL 



/ Kn{x, 



x) dx + ( I Kn[x,x)dx 



J JliXL 



y) dxdy 



It now follows that 



Var(#/,) = f Kn{x, x) dx - /"/" kI{x, y) dxdy. 

To get a more convenient formula to work with one can now use the identities 
Kn{x,y) = Kniy,x) and 



/ 



Kn{x, y)Kn{y, z) dy = K{x, z) 



to get 



Var(#/,)= / (/ Kl{x,y)dY\ dx - // Kl{x,y) dxdy 

Jli \JR J J Jli-Kli 

= 11 Kl{x,y) dxdy. 

JJliXlf 

In more generality one gets 

m 

E[ai#Ii + . . . + am#Im] = yZ(^i Kn{x, x) dx 
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and 

(ai#/l + . . . + am#Ira? 



m /I 



j=l \A;=1 / i^j \k=l J \k=l J 

=: Si + S2. 

Prom the calculations above we know that 
E[Si] = X { / ^n{x, x) dx + (^J^ Kn{x, x) dx^ 

-// Klix,y)dxdy 

J JliXli 

so it remains to calculate E[52]. We have 

(n \ / ^ \ 

^Xii{xk)\ [^Xij{xk)\ =^Xii{xk)Xh{xi) 
k=l ) \k=\ J ki-l 



and hence 



^[-5*2] = X / / Pn,2{x, y) dxdy 
= y2aiaj[ / Kn{x,x)dx Kn{x,x)dx- // i^^(x,y)dxdy . 



^+3 

Since 

(E[ai#7i + . . . + a„#7^])' 

2 m 



= V^ciM / ii'„(x,x)dx) + OiOj / i^n(a^, a;) dx / Er„(x,x)dx 



we finally get (with manipulations as before) 
Var(Qi#/i + . . . + am#Im) 

= ^^i ^n{x, y) dxdy - ^ aiaj II Kl{x,y) dxdy. 
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